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ABSTRACT 

The observations of a nearby low-luminosity gamma-ray burst (GRB) 060218 associ- 
ated with supernova SN 2006aj may imply an interesting astronomical picture where 
a supernova shock breakout locates behind a relativistic GRB jet. Based on this pic- 
ture, we study neutrino emission for early afterglows of GRB 060218-like GRBs, where 
neutrinos are expected to be produced from photopion interactions in a GRB blast 
wave that propagates into a dense wind. Relativistic protons for the interactions are 
accelerated by an external shock, while target photons are basically provided by the 
incoming thermal emission from the shock breakout and its inverse-Compton scat- 
tered component. Because of a high estimated event rate of low-luminosity GRBs. 
we would have more opportunities to detect afterglow neutrinos from a single nearby 
GRB event of this type by IceCube. Such a possible detection could provide evidence 
for the picture described above. 
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1 INTRODUCTION 

Gamma-ray burst (GRB) 060218 associated with SN 2006aj 
discovered by Swift (Campana et al. 2006) provides a new 
example of low- luminosity GRBs (LL-GRBs), as its isotropic 
equivalent energy (~ 6 x 10*^ erg) is 100 to 1000 times less 
but its duration (Tgo = 2100 ± 100 s) is much longer than 
those of conventional high- luminosity GRBs. More interest- 
ingly, besides an usual non-thermal component in its early 
X-ray spectrum, a surprising thermal component was ob- 
served by the Swift XRT during both burst and afterglow 
phases. Fitting with a blackbody spectrum, the temperature 
of this thermal component was inferred to be kT ~ 0.17 keV 
during the first 3 ks. When t >10 ks, however, the peak en- 
ergy of the blackbody decreased and then passed through 
the Swift UVOT energy range at ~100 ks (Campana et al. 
2006; Blustin 2007). 

To explain the prompt emission, in principle, a model 
based on the internal dissipation of relativistic ejecta may 
be valid in the case of GRB 060218. The relativistic GRB 
ejecta, which interacts with a dense wind surrounding the 
progenitor, has also been required to understand a power- 
law decaying afterglow in X-ray and radio bands about ~10 
ks after the burst, although some complications beyond the 
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standard afterglow model are involved (Soderberg et al. 
2006; Fan et al. 2006; Waxman et al. 2007). Furthermore, as 
suggested by Wang & Meszaros (2006) , the soft X-ray ther- 
mal emission could arise from a shock breakout, namely, a 
hot cocoon that breaks out from the supernova ejecta and 
the stellar wind. In detail, a more rapid part of a jet mov- 
ing in the envelope and the dense wind is accelerated to 
a highly-relativistic velocity to produce the GRB, while a 
slower part of the jet together with the outermost parts of 
the envelope becomes a mildly relativistic cocoon (Meszaros 
& Rees 2001; Ramirez-Ruiz et al. 2002; Zhang et al. 2003), 
which locates behind the GRB blast wave. 

Although the GRB blast wave runs in front of the shock 
breakout, the thermal emission from the latter outshines 
the former persistently until the emission of the breakout 
is switched off. Thus, the emission properties of the GRB 
blast wave (consisting of external shock-accelerated elec- 
trons and protons) should be influenced by the incoming 
thermal photons signiflcantly during both burst and early 
afterglow phases. On one hand, the cooling of the relativis- 
tic electrons, which upscatter the thermal photons, could be 
dominated by inverse-Compton (IC) radiation rather than 
synchrotron radiation (Wang & Meszaros 2006). On the 
other hand, inferred from the observations, the intensity of 
the thermal emission could be comparable in the same band 
to the one of the prompt emission due to internal dissipa- 
tions and much larger than the one of the afterglow emission 
due to an external shock. Therefore, the thermal photons as 
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target photons for photopfon interactions could also play an 
important role in the energy loss of the relativistic protons 
and thus influence or even dominate neutrino emission of 
the GRB blast wave. 

It has been widely studied that conventional GRBs in 
the standard internal-external shock model emit high energy 
neutrinos during the burst, early afterglow, and X-ray flare 
phases (Waxman & Bahcall 1997, 2000; Dai & Lu 2001; Der- 
mer 2002; Dermer & Atoyan 2003; Asano 2005; Murase & 
Nagataki 2006a, 2006b; Murase 2007; Gupta & Zhang 2007). 
In contrast to the conventional GRBs, LL-GRBs may have a 
much higher event rate (several hundred to thousand events 
Gpc~^yr~^), which is inferred from the fact that two typ- 
ical nearby LL-GRBs, i.e., GRB 060218 and GRB 980425, 
have been observed within a relatively short period of time 
(Cobb et al. 2006; Soderberg 2006; Liang et al. 2007). The 
high event rate implies that the contribution from LL-GRBs 
to the diffuse neutrino background may be important and 
that we have more opportunities to detect neutrinos from 
a very nearby single LL-GRB event. Therefore, Murase et 
al. (2006) and Gupta & Zhang (2006) recently studied the 
neutrino emission properties of LL-GRBs during their burst 
phase using the internal shock model, in which the target 
photons for photopion interactions are mainly provided by 
internal shock-driven non-thermal emission. In this paper, 
however, we will focus on the early afterglow neutrino emis- 
sion. During this phase, relativistic protons are accelerated 
by an external shock (rather than internal shocks) and tar- 
get photons are dominated by the incoming thermal emis- 
sion and even its IC scattered component. 

This paper is organized as follows. In section 2 we briefly 
describe the dynamics of a GRB blast wave propagating into 
a surrounding dense wind. In section 3, we give the photon 
distribution in the blast wave by considering the incoming 
thermal emission and its IC scattered component, but the 
weak synchrotron radiation of the electrons is ignored. In 
section 4, neutrino spectra are derived formally with an en- 
ergy loss timescale of protons due to photopion interactions. 
In section 5, we calculate the timescale and then the neu- 
trino spectra using the target photon spectra obtained in 
section 3 and an experiential fltting formula of the cross 
section of photopion interactions. In addition, the peak of 
a neutrino spectrum is also estimated analytically by using 
A— approximation. Finally, a summary is given in section 6. 



2 DYNAMICS OF A GRB BLAST WAVE 

We consider a GRB jet with isotropic equivalent energy 
E = 10^°-B50 erg (hereafter Q^, = Q/10^) expanding 
into a dense wind medium with density proflle p(r) = 
Ar~^. Here, the coefflcient A is determined by the mass 
loss rate and velocity of the wind of the progenitor, i.e., 
A = M/4tvv^ = 5.0 X 10"gcm"M,, where A, = 
[A//(10"^Mq yr-^)]K/(10^km s-^)]-\ From Dai & Lu 
(1998) and Chevalier & Li (2000), we get the Lorentz factor 
and radius of the GRB blast wave (i.e., the external-shocked 
wind gas) respectively as 
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They satisfy r = ST^ct, which gives rise to a relationship, 
t' = (16/3)rt, between the dynamic time t' measured in the 
rest frame of the blast wave and the observed time t (Dai & 
Lu 1998). 

As the circum-burst wind materials are swept up and 
shocked, most of the heated electrons before cooling con- 
centrate at the minimum Lorentz factor 7e.m ~ ^e^F — 

659 ee,-iEl^''A-^^%^^^. The symbol = e,{p-2)/{p~l), 
where ee is the usual equipartition factor of the hot electrons 
and p is the electron's energy distribution index (where p > 2 
is only considered). Meanwhile, a fraction es of the inter- 
nal energy is assumed to be occupied by a magnetic field, 
and then the strengthen of the magnetic field is calculated 
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by B' ~ (327resr>c^)^''^ = 78G e'J'_^E, 
nally, the other energy (a fraction of tp = 1 — te — eb) is car- 
ried by the accelerated protons. For these protons, we can 
estimate their maximum energy by -Ep.max = 2eB'r/3 = 
4.8 X 10^°GeV e]^^_iEll,*AV%^^'^ by equating the accel- 
eration time to the shorter of the dynamic time and the 
synchrotron cooling time (Razzaque el al. 2006). However, 
the minimum energy of the protons is unknown, but the 
corresponding Lorentz factor 7p_niin is thought to be close 

to ~ r. 



3 PHOTON EMISSION 

The photons in the GRB blast wave have two origins, i.e., the 
blast wave self and the inner supernova shock breakout. The 
electrons in the blast wave emit photons via synchrotron and 
IC scattering processes. Moreover, as analyzed by Wang & 
Meszaros (2006), the synchrotron radiation (peaking within 
X-ray band) of the blast wave electrons is inferred from the 
observations to be much weaker than the incoming thermal 
emission, and thus the cooling of the electrons should be 
dominated by their IC scattering off the thermal photons. 
Therefore, in following calculations, we consider the thermal 
emission and its subsequent IC scattered component only. 

The properties of the supernova shock breakout have 
been unclear to date. We suppose that it has a constant 
blackbody temperature of kT = 0.1keV(fcr)_i and a con- 
stant radius of the emission region of J? = 10^^cm7?i2. 
The lifetime isB of this high-temperature emission is about 
thousands of seconds, which is considered to be several 
to several ten times longer than the duration of the 
GRB. Then, the isotropic equivalent luminosity and en- 
ergy of the shock breakout can be estimated by Lsb ~ 
A-kR'cjT'^ = 1.3 X lO^^erg s"^ {kT)t^Rl2 and _Bsb = 1.3 x 
10'**erg (fcT)li_Ri2isB,3, respectively. Meanwhile, it is easy 
to write the monochromatic number density of these thermal 
photons at the breakout as 
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where the function 4>{x) — x^/{e^ ~ !)• Assuming that the 
photons propagate freely before they reach the GRB blast 
wave at radius r, we can calculate the density of the thermal 
photons in the blast wave by multiplying a factor [R/r)'^ 
to Eq. ((Sjl. Subsequently, after Lorentz transformation, we 
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obtain the density of the incoming photons in the blast wave 
measured in its rest frame by 
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where i?^,,pki = 3fcT/r is the peak energy of the black body 
spectrum. When these photons cross the blast wave, a part 
of them should be upscattered by the relativistic electrons. 
The energy of the IC scattered photons can be estimated by 
E'ry ic = 27'^ mE'h-i s-nd the corresponding density by 
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where £^7,pk2 = 67'e ^kT/F. The probability of the scatter- 
ing is represented by the photon optical depth of the blast 
wave, T = aT{A/mpr) = 6.5 x lO"'^ E^o^'^Al^^t^^^^, where 
ITT is the Thomson cross section. According to this esti- 
mation, Wang & Meszaros (2006) predicted that the early 
afterglow spectra of GRB 060218-like GRBs may have a bi- 
modal profile peaking at 



S^.pki = 0.3 keV (fcT)_i 
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Thus, a significant sub-GcV or GeV emission component 
accompanying the thermal emission would be detectable 
with the upcoming Gamma-ray Large Area Space Telescope, 
which could provide evidence for the GRB jet. 



4 NEUTRINO PRODUCTION 

Since relativistic protons in the GRB blast wave are im- 
mersed in the photon field described above, the protons 
would lose their energy to produce mesons such as tt*^ and 
TT* etc, and subsequently generate neutrinos by the decay of 
71*, i.e., 71* + Uf_,{v^) e"^ + Ue{i'e) + t^^ + iyfj,. During 

these processes, the energy loss rate of a proton with energy 
E'p ~ "ipmpC? can be calculated by (Waxman & Bahcall 
19970 
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where is the cross section of photopion interactions 

for a target photon with energy E in the proton's rest frame. 



^ To obtain this expression, an isotropic target photon field is 
required. However, in our model, the radially incoming photon 
field is seen by the protons in the blast wave anisotropically. This 
gives an extra complication for a more realistic consideration. For 
simplicity, we ignore the anisotropic effect in our calculations. 



^ is the inelasticity defined as the fraction of energy loss of 
a proton to the resultant pious, and Eth = 0.15GeV is the 
threshold energy of the interactions. Equation ^ yields that 

the energy of the protons decreases as exp — [dt' /t' -k) 
In our scenario, if the shock-breakout emission could last for 
a period of tsB, the fraction of the energy loss of the protons 
to pions could be calculated by 



/jT = 1 - exp - / — 
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where fgg = (16/3)rtsB. In order to calculate t'^, the crucial 
input in the model is the target photon spectrum n'{E'j). 
From Eqs. Q and we know that n'{E!y) depends on 
both r and F and thus the value of t'^ could evolve with 
time. However, if t'^ is independent of time or varies with 
time slowly, Eq. Q can be also approximated by 
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as usual, especially for analytical calculations. 

To be specific, the energy loss of the protons is shared by 
71* and 71" with a certain ratio. Unfortunately, it is not easy 
to fix this ratio due to the complications arising from various 
single-pion and multipion production processes. In following 
calculations, we simply take it to be a constant, 71* : 71" — 
2 : 1, as in Asano (2005). Furthermore, two resultant muon- 
neutrinos from the decay of a 71* could inherit half of the 
pion's energy roughly evenly. Therefore, we can relate the 
neutrino energy Ei, to the energy loss of the primary proton 
by 
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and give an observed time-integrated muon-neutrino spec- 
trum by 
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where Di is the luminosity distance of the burst. As usual, 
we assume the energy distribution of the shock-accelerated 
protons to be (dNp/dEp) oc Ep^, where the proportional 
coefficient can be calculated by tpi?/ ln(i?p,niax/i?p,miii). 

In addition, because of the presence of the magnetic 
field, the ultrahigh energy pions and muons would lose their 
energy via synchrotron radiation before decay. This leads to 
breaks in the neutrino spectrum at (Murase 2007) 
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where t^, = 2.6 x 10~**s and = 2.2 x 10~''s are the 
mean lifetimes of pions and muons in their rest frames. 
Above El I, , the neutrino flux would be suppressed by a 
factor {El, / E1_^)~'^ (Rachen & Meszaros 1998; Razzaque et 
al. 2006). However, as pointed out by Asano & Nagataki 
(2006), neutral kaons can survive in the magnetic field, while 
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the ultrahigh-energy charged pions and muons cool rapidly. 
Moreover, because kaons have a larger rest mass than pions 
and muons, charged kaons can reach higher energy although 
they also suffer from synchrotron cooling. Thus, decay of 
kaons, which is not taken into account in our calculations, 
may dominate neutrino emission above ~ 10* — lO^GeV. 

Now, by inserting Eqs. ([4]) and ([U into Eq. ([U and 
then into Eq. (|9]) to get we can easily obtain the ob- 
served neutrino spectra from Eq. (|12|l for our scenario. The 
remaining task is only to express the cross section cr,r(£^) 
and inelasticity for photopion interactions. 



5 RESULTS 

Since the cross section of photopion interactions peaks at 
i?A — 0.3GeV due to the A(1232)-resonance, the integration 
over E in Eq. ((Sjl can be roughly approximated by 

f oc 

t'-^ ~ T^(^.,aUE^5E / n'{E'^)E':,^dE'-^, (15) 

where cr-^,A ~ 0.5mbarn, ~ 0.2, and the peak width is 
about SE « 0.2GeV. Inserting Eq. (g]) or (0) into Eq. ([TSll . 
we use the approximative formula = min[tgB/t^,l] to 
obtain 

.16 R^8nk^2E:^ - 
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where ? = 1 and T/{2'y' ^ ^) for pre- and post- 
upscattered target photons, respectively. The dimension- 
less variable e* is defined by £« = 3E/\/{2"f'pE'^ py.) = 
3£,A^^ EAmpC^ / (SE^E-y^pii). In the case of frr < 1, the peak 
value of /,r reading 

A.pk - 3?— j^a^M^SETtsB (17) 

is at e* = 1.8, which gives rise to the relationship between 
the peak energies of the neutrino and photon spectra as 
£'i/,pk£'7,pk = O.Oir^GeV^. Considering the bimodal distri- 
bution of the target photons peaking at -E^-^pki and -E-y,pk2, 
two peaks are also expected in the resultant neutrino spec- 
trum but only the one determined by i?-y,pki could fall into 
the high energy range {E^ > TeV) of our interest at 
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In other words, the target photons for photopion interactions 
of interest are contributed by the incoming thermal emission 
mainly. The value of the differential neutrino fluence at Ev,pk 
reads 

[S,^(?i,]pk = 2.0 X lO-'^erg cm"' ep(fcr)lii??2^. AT25.5- (19) 

which is calculated by using the peak value of as 

A,pk = 0.02 {kTyURl^E^o^A,. (20) 

On the other hand, when /,r = 1, E^cp,^ would reach an 
upper limit as 1.2 x 10"'' erg cm~^ epE^oD^^s 5) which is 
determined by the total energy carried by the protons in 
the GRB blast wave. 



Although it is convenient and effective to use the 
A— approximation to estimate the peak of a neutrino spec- 
trum, the A— approximation would lead to an remarkable 
underestimation of the neutrino flux above the peak energy 
due to the non-zero cross section of photopion interactions 
in high energy regions. So, for more careful calculations, we 
provide an experiential fitting formula for the cross section 
as shown in Eq. (|A5|) . which is extrapolated from experimen- 
tal data taken from particle data group (Yao et al. 2006). 
However, since we can not find a simple expression for the 
inelasticity, we take ^ = 0.2 for all energy regions roughly, 
which may leads to a mild underestimation of the neutrino 
flux in the high energy regions. Finally, with these inputs, 
we plot the observed time-integrated muon-neutrino spectra 
in Fig. 1. Obviously, two plateaus exist in the neutrino spec- 
tra. To be specific, as shown in the upper panel of Fig.l, the 
high-energy plateau is produced by the lower energy ther- 
mal photons, while the low-energy plateau is produced by 
the higher energy IC scattered photons. In addition, from a 
comparison shown in the lower panel of Fig. 1, we can see 
that the approximation for /,r in Eq. pop is feasible to some 
extent for the thermal seed photon-dominated photopion in- 
teractions, but not for the IC scattered photon-dominated 
interactions. This difference of these two kinds of interaction 
arises from different temporal behaviors of t'^. 

Next let's discuss the detectability of the afterglow neu- 
trinos, using the following fitting formula for the probability 
of detecting muon-neutrinos by IceCube (loka et al. 2005; 
Razzaque et al. 2004) 
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where /3 = 1.35 for < 3.2 x 10*GeV, while /3 = 0.55 
for E,y > 3.2 X lO^GeV. The number of muon events from 
muon-neutrinos above TeV energy is given by 



t>uPvdEv, 



(22) 



where Adct ~ Ikm^ is the geometrical detector area. In- 
serting Eqs. (|12|) and (|21|l into the above integral, we ob- 
tain Nf_i ~ 0.1 for the parameter set (-E50 = 1, ^4* = 10, 
(fcT)_i = 2, i?i2 = 1 and tsB.a = 3) inferred from GRB 
060218 for a very nearby LL-GRB at 50 Mpc, where a LL- 
GRB event is expected to be observed within a many-years 
observation. According to this estimation, we expect opti- 
mistically that IceCube may be able to detect afterglow neu- 
trinos from one LL-GRB event in the following decades. If 
such a detection comes true, the afterglow neutrino emission 
accompanying the soft X-ray thermal and sub-GeV or GeV 
emissions from a GRB 060218-like GRB event would pro- 
vide strong evidence for the picture that a supernova shock 
breakout locates behind a relativistic GRB jet, and further 
would be used to constrain the model parameters severely. 

Besides the possible detection of neutrinos from a single 
LL-GRB event, the contribution to the neutrino background 
from LL-GRBs is also expected to be important. We can es- 
timate the diffuse muon-neutrino fiux arising from afterglow 
neutrino emission of LL-GRBs by (Waxman & Bahcall 1998; 
Murase et al. 2006) 
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Figure 1. The time-integrated afterglow muon-neutrino {v/^+Pi^) 
spectra for one GRB event. The sohd hnes are calculated by using 
the expressions for in Eq. ||9} and for cr^{E) in Eq. IIA5I I. 
Upper panel: The contributions to the total neutrino emission by 
the two target photon components are represented by the dashed 
and dotted lines, respectively. Lower panel: The dashed line is 
obtained by an approximation for the time integration as in Eq. 
moi l, and the peak estimated by the A— approximation is labeled 
by an open circle. In all cases, we take the model parameters 
E50, A«, (fcT)_i, -R12, Di 25.5 and <sb 3 to be unity and the 
equipartition factors Ce = 0.3, es = 0.1 and thus ep = 0.6. 



= 2.5 X 10""GeV cm"^s"^sr"^ ep(fcr)lii??2^ 
/ i?LL(0) \ ff, 



V 500Gpc-»yr-i 
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where Hq — 71km s^^Mpc^^, is the beaming factor, and 
is the correction factor for the possible contribution from 
high-redshift sources. In the above estimation, the approxi- 
mative value of in Eq. (|20|l is applied. By comparing Eq. 
(f23)) to Eq. (3) of Murase et al. (2006), we find that, for LL- 
GRBs, the contribution to the diffuse neutrino background 
by the early afterglow neutrino emission may be relatively 
smaller than or even comparable to (e.g., for model param- 
eters ep = 0.6, {kT)-i = 2, R12 = 1, and A, = 10) that by 
the burst neutrino emission. 

Finally, we would like to refer the reader to neutrino 
oscillation, which will change neutrino flavor ratio from Ve '■ 
: ^'t- — 1 : 2 : at the source to 1 : 1 : 1 at the earth. 
This thus leads to the fact that the observed muon-neutrino 
fluxes estimated above should be reduced further by a factor 
of ~ 2. 



6 SUMMARY 

The surprising soft X-ray thermal emission during both 
burst and afterglow phases of GRB 060218/SN 2006aj was 
proposed to be due to the breakout from a strong stellar 
wind of a radiation-dominated shock. This shock breakout 
was further thought to locate behind a relativistic GRB 



jet, which is required by understanding the burst emission 
and the power-law decaying afterglow emission. Wang & 
Meszaros (2006) suggested that a sub-GeV or GeV emission 
produced by IC scattering of the thermal photons by the 
relativistic electrons in the GRB blast wave could give evi- 
dence for this astronomical picture. In this paper, we stud- 
ied another possible implication, namely, afterglow neutrino 
emission. The neutrinos are produced by photopion interac- 
tions of relativistic protons, which could be accelerated by a 
relativistic external shock. The target photons in the interac- 
tions are contributed by the incoming thermal emission and 
its upscattered component. By considering the high event 
rate of LL-GRBs, we argue optimistically that the afterglow 
neutrinos from very nearby (several tens of Mpc) LL-GRBs 
may be detected by IceCube in the following decades. We 
believe the detection of these expected afterglow neutrinos 
is helpful to uncover the nature of GRB 060218-like GRBs. 
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APPENDIX A: CROSS SECTION FITS 

The flts to the total cross section for P7 interactions have 
been widely studied (e.g., Rachen 1996; Miicke et al. 2000). 
In physics, this cross section is contributed by resonant ex- 
citations and direct (non-resonant) single-pion production 
processes in the resonant energy regions and by statistical 
multipion production processes mainly and diffractive scat- 
tering slightly in the high energy region [Ea >0.727 GeV). 

It is known that the cross section for a resonance is 
given by the Breit-Wigner formula (Miicke et al. 2000) 



aoF^s 



£2 (s-mlc-^y + r^s 



(Al) 



where ^/s = {rUpC^ + 2mpc'^ E)^^^ is the total energy of the 
colliding photon and proton in the mass-center frame, and 
m_R and f are the nominal mass and width of the reso- 
nance, respectively. The coefBcient tro is determined by the 
resonance angular momentum and the electromagnetic exci- 
tation strength. For nine important resonances in p7 inter- 
actions, we take the related parameters from Miicke et al. 
(2000) and list them in Table 1. Then, the total cross section 
contributed by these resonances can be written as 



1 — exp 



0.15 



0.08 



ME), 



(A2) 



where the suppression factor in the square bracket repre- 
sents the threshold i5th = 0.15GeV of P7 interactions with 
X = E/GeV. Moreover, Rachen (1996) found that the cross 
section in the high energy region {Eq >0.727 GeV) can be 
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Table Al. Parameters for resonances. 



Name 


mflcVGeV 


r/GeV 


(TO / /ibarn 


A(1232) 


1.231 


0.11 


31.125 


Af(1440) 


1.44 


0.35 


1.389 


Af(1520) 


1.515 


0.11 


25.567 


Ar(1535) 


1.525 


0.1 


6.948 


Af(1650) 


1.675 


0.16 


2.779 


Af(1680) 


1.68 


0.125 


17.508 


A(1700) 


1.69 


0.29 


11.116 


A(1905) 


1.895 


0.35 


1.667 


A(1950) 


1.95 


0.3 


11.116 



fitted bj0 



1 — exp 



0.727- 



0.8 



(0.067y""*^ +Q.12^V 



X 

0.453\ 



mbarn. 



(A3) 



where y — s/GeV^. Subtracting tlie two components ex- 
pressed by Eas. lA2l and lA3l from the experimental data taken 
from particle data group (Yao et al. 2006), we find the resid- 
uals to the total cross section exhibit a broken power-low 
behavior, which yields 

'0.15 

exp ( 



<T3(S) 



0.072a:" 



0.08 



+ 1 



+ 1 



+ 1 



(A4) 



mbarn. 



Roughly speaking, the above formula could be related with 
the direct single-pion production processes. Finally, com- 
bining the three components, we can express the total cross 
section of P7 interactions by 



tJi + as, for 0.15GeV < Eq < 0.727GeV; 
ai+a2+ as, for Eq > 0.727GeV. 



(A5) 



We confront this experiential formula with the experimental 
data in Fig. Al. It can be seen the fits is good, although the 
formula can not describe the detailed physics. 
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